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I: 
o 
in 

A simple lattice model describing the recombination dynamics in visible 
light emitting porous Silicon is presented. In the model, each occupied lattice 
site represents a Si crystal of nanometer size. The disordered structure of 
porous Silicon is modeled by modified random percolation networks in two 
5_i ■ and three dimensions. Both correlated (excitons) and uncorrelated electron- 

hole pairs have been studied. Radiative and non-radiative processes as well as 
hopping between nearest neighbor occupied sites are taken into account. By 
means of extensive Monte-Carlo simulations, we show that the recombination 
dynamics in porous Silicon is due to a dispersive diffusion of excitons in a 
disordered arrangement of interconnected Si quantum dots. The simulated 
luminescence decay for the excitons shows a stretched exponential lineshape 
while for uncorrelated electron-hole pairs a power law decay is suggested. Our 
results successfully account for the recombination dynamics recently observed 
in the experiments. The present model is a prototype for a larger class of 
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models describing diffusion of particles in a complex disordered system. 
Pacs: 71.55.Jv;72.20.jv;78.47.+p;61.43.Bn 
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I. INTRODUCTION 



Silicon, the most studied semiconductor, is not a good light emitter, especially in the 
visible range, because of its indirect band gap transition in the infrared (1.1 eV). The re- 
cently discovered optical properties of porous forms of Silicon (p-Si) have therefore attracted 
considerable interest. [l|-|4j In addition to the internal nanometer sized structures which are 
distributed in space according to a complex topology, || these new materials have similar 
spectral behaviors, the most prominent are a rather high luminescence quantum efficiency 
and an unexpected wide range of luminescence lifetimes (ranging from microseconds to mil- 
liseconds), depending on temperature. 

The mechanism of the luminescence emission is still unclear. Three main models have 
been proposed and survived to the recent discussions: 

1. quantum recombination model, 

2. surface state model, M 



3. molecular recombination model. 10-12 



The first two models agree on the fact that the energy spectrum of p-Si is a result of 
the quantum confinement of carriers in Si nano-crystals. Still under debate, is the shape 
of the nano-crystal that has been described either in terms of "undulating" quantum wires 
or quantum dots. || The quantum recombination and the surface state models differ in 
their predictions about the origin of the luminescence. In the first case, the luminescence 
is due to confined and localized excitons, while in the second, the luminescence is due to 
electron-hole pair recombination, where the electron and the hole could be found either in a 
bulk nano-crystal state ('extended' state), or being trapped into a surface nano-crystal state 
('localized' state). The third model, the molecular recombination model, considers that some 
particular molecular species like polysilane chains (e. g. SiH chain [|10j) or siloxene-like rings 
(Si-O-H) JTTJ] form on the surface of the Si nano-crystals. The luminescence is due to the 
carrier trapping into these species. 



Recent experiments [|13| -|i8|,?,|20|| have provided clear evidence of anomalous relaxation 



behavior of the luminescence. The decay lineshape, for a single observation energy, is not 
described by a single exponential function. This non-exponential behavior can be commonly 
described by a stretched exponential function, |pT|-p3| defined as 



L(t) = L(0) exp [- (t/rf] (1) 

where L(t) is the time dependent luminescence intensity, r is a lifetime and (3 < 1 a dispersion 
exponent. In general, values of /3 < 1 correspond to the existence of a broad distribution of 
lifetimes. In some circumstances, such broad distributions may be the result of a diffusive 
motion of the excited carriers. This may be the case in p-Si, as suggested recently in terms 
of trap-controlled hopping processes. |18| We elaborate this idea further in this work. 



Motivated by these experimental results, and by the fact that little is known about the 
recombination dynamics of the charge carriers in these materials, we have initiated a detailed 
numerical study of the underlying transport behavior in p-Si by means of Monte-Carlo (MC) 
simulations. 

In this work, we present a simple model of p-Si in which nanometer sized crystals (nano- 
crystals), characterized by a distribution of radiative and non-radiative recombination times, 
are assumed to be randomly placed at the sites of percolation-like clusters defined on square 
and simple cubic lattices. Charge carriers are allowed to hop between nearest-neighbor 
occupied cluster sites. The competing effect between radiative and non-radiative transitions 
in a single Si nano-crystal, as well as the effects of geometrical constraints on hopping of 
carriers due to the complex topology of percolation clusters, are discussed. 

Even though our simulations are based on the quantum recombination model, |||7| they 
are consistent with other possible recombination mechanisms because the model concentrates 
on the dynamics of the recombination and not on the recombination process itself. However, 
we found that most of the data are better explained by our simulations when the quantum 
recombination model is assumed. In addition, since our model is based on quite general 
assumptions about the geometrical disorder involved, we expect that it may be applied to 



a larger class of physical systems characterized by both a structural and local disorder, of 
which p-Si is just one possible example. 

The paper is organized as follows: in Section II the model is described, and exact so- 
lutions obtained in the case of isolated Si nano-crystals are discussed. In Sections III and 
IV, simulation results for one-, and two- and three-dimensional systems, respectively, are 
presented. In Section V, a discussion of the present results is presented and compared to 
other results already known in the literature. Finally, in Section VI we give some concluding 
remarks. 

II. MODEL FOR A SI NANOMETER SIZED CRYSTALS 

We consider a lattice model for p-Si in which occupied and empty sites are present. 
Each occupied site of the lattice represents a Si nano-crystal. In one dimension, all the 
sites are occupied by Si nano-crystals, while in higher dimensions the occupied sites are 
randomly distributed in space and interconnected to each other according to the topology 
of percolation-like clusters (see Section IV). We further assume that the physical properties, 
relevant for our problem, of such Si nano-crystals are suitably described in terms of quantum 
dots (QD), for which many theoretical results are presently available. 

A. The size and radiative energy of a Si nano-crystal 

Each Si nano-crystal, or equivalently, each occupied site of the lattice, is characterized 
by a size d, representing its mean diameter (typically in the range 15 to 55 A p4]j). Since 
the actual spatial distribution of nano-crystal sizes is not known a priori, the sizes di at 
site i are chosen randomly according to a Gaussian distribution P(d). We used P(d) ~ 
(1/od) exp[— (l/2<7^) (d — d ) 2 ], centered at do = 32.4 A and having a width ad = 0.8 A. 

In one dimension, this local randomness is the only source of disorder in the model, while 
in higher dimensions, an additional (geometrical) disorder, due to the complex structure of 
percolation-like clusters, is also present (see Section IV). 
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Due to the experimental method of preparation of p-Si, however, some correlation be- 
tween the sizes d may actually be present for nearby Si nano-crystals, from which one may 
reasonably expect that strong local variations in d are not likely to occur. || To take this 
effect into account in a simple way, the sizes d corresponding to nearest neighbor occupied 
sites are not allowed to differ by more than a prefixed amount. We are going to specify this 
point below. 

By neglecting the excitonic binding energy, the optical band-gap energy E associated 
to a given Si nano-crystal is assumed to correspond to the emission energy determined 
experimentally. Ei is related to the size di through the power-law relation recently suggested 
by calculations performed for spherical quantum dots. [26| Based on these calculations, we 



assume Ei = E (d /di) n , with n = 1.39, where E = 1.86 eV is the transition energy for a 
Si nano-crystal of size do. The corresponding distribution for E is then of the form, 

P{E) ~ ± exp[-(dS/20*) ((E /E)^ - I) 2 ] (2) 

centered around Eq and having a width at half-maximum of about 0.16 eV (see Fig. [I]). 

It is instructive to obtain estimates of the total density of states p(e) resulting from 
these distributions. For simplicity, an ensemble of quantum dots of cubic shape and side 
L x = L y = L z = d is considered. |27| The (confined) electronic states have energies e(n X)J/)2 ) = 
£o(nx + n y + n z ), where eo = h 2 n/ (2/j,d 2 ), n x , n y and n z are the principal quantum numbers 
for motion in the x, y and z directions, n X)V>z > 1 and fi is the effective mass. The transition 
energy E can be related to the ground state energy of the QD by E = 3e + -EWik, where 
-^Buik = 115 eV is the bulk-Si indirect gap. The density of states for a single quantum dot 
of size d with an energy barrier of height V(QD) is then |27j 



PM - £ 5(e - e(n x , y , z ) - E Bnlk ), (3) 



where e is an energy and the total density of states is 



p(c)=E P (d) Pd (e). (4) 
d 
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This quantity is plotted in Fig. [I] assuming a Gaussian broadening of the single QD state 
of 0.01 eV. The density-of-states of an ordered ensemble of QDs should reveal a series 
of distinct peaks corresponding to the different confined states e{n x ,y,z) in the QD. The 
random Gaussian distribution of sizes, i. e. E, smears out the discrete peaks and yields 
an approximately exponential behavior of the form p(e) ~ exp(e/e a ), with e a = 0.35 eV, at 
large e. 

To deal with the above mentioned spatial correlations in d, the difference in energy 
between nearest-neighbor sites i and j, A£?y = Ej — E iy is allowed to take absolute values 
smaller than a given cutoff E cnt . This is the criterion employed to accept a given value of dj, 
with respect to the previous value at site i, from the Gaussian distribution. In dimensions 
higher than one, this criterion is applied simultaneously with the growing process (see Section 
IV) . For convenience, we consider variations in energy, rather than variations in size to select 
the values of the local quantities E and d. The results are not sensitive to this choice. Typical 
calculations have been performed for values in the range 10 -4 < E cut < 0.4 eV (see Section 
III). E cut is the first parameter in the model. 



Within each Si nano-crystal, both radiative and non-radiative recombination processes 
occur when excited electron- hole pairs are present. Both processes depend on the size d of 
the nano-crystal. 

The emission process is assumed to be due to excitonic recombinations in a QD. In this 
case, it was shown that the radiative recombination time, r ra d, results from the thermal 
balance between the occupation of the exchange splitted singlet and triplet excitonic states, 



T S i n g (the radiative lifetime for the singlet excitonic state) is obtained from the values calcu- 



B. The radiative and non-radiative recombination times 





lated in Ref. || 
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where n s = 1.5, and r^ g is the second parameter in the model. At the present, no theoretical 
calculations are available for the radiative lifetime of the triplet excitonic state, r trip i. Thus, 
Ttripi is estimated from the experimental data reported in Ref. [[7], 

T tri pl(£) = r t ( r t [T sing (E)f° (7) 

where r^ pl = 2368 /zs, and a s = 0.307. 

Finally, the quantum confinement enhanced exchange energy AE X = £f {djdf, where 
Ef = 3.165 x 10" 4 eV is the exchange energy in crystalline Si and d x = 43 A, is obtained 



from the calculations in Ref. 29 



Non-radiative recombinations have been considered as multiphonon transition processes 



for trapping on a deep center, J3(J for which the non-radiative recombination time, r nr , is 
given by 

p 2n T 5 

where tit = [exp(hu ph J k B T) — 1] _1 is the density of phonons with energy hu ph , and S and p 
are related to the details of the capturing center. Here we have used nano-crystal size inde- 
pendent quantities, with huph/ks = 800 K, S = 1, and p = 25. If bulk-Si optical phonons 
are involved, Tiuj ph /k B ~ 1100 K; lower values can be justified by the typical softening of the 
phonon modes in a confined system on a localized center. t$ is the third parameter in the 
model. 

It should be emphasized that the actual nature and origin of non-radiative phonon tran- 
sition processes in p-Si nano-crystals remain to be understood. Here, such processes are 
included because they certainly play a very important role in the recombination dynamics 
of carriers. [§,0 The form for r nr in Eq. |8| can be only considered as a first approximation to 
the actual form of the non-radiative recombination time, yet it is expected to describe the 
physics involved correctly, at least in a qualitative fashion. In this work, the same r nr has 
been employed for all the Si nano-crystals. 
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C. The hopping rates between Si nano-crystals 



In addition to the internal structure of Si nano-crystals, hopping of electrons, holes 
and excitons, may actually occur between nearby nano-crystals. To describe the diffusion of 
carriers in p-Si, we consider hopping processes only between nearest neighbor (n.n.) occupied 
sites of the lattice for two different models. In the first model, only one type of carrier 
hops, representing the exciton, while in the second, two different types of carriers can hop, 
representing the electron and hole. If the electron and hole (e-h) motions are strongly 
correlated, e.g. if the hole closely follows the electron as it moves in the system, one has 
essentially the first model. As we will see in Sections III and IV, a qualitatively different 
behavior is obtained when e-h hops are completely uncorrelated. In both models, correlated 
and uncorrelated e-h pairs, the on site radiative recombination is assumed to be excitonic. 

The transition rates for a carrier of type x (x = ex, representing the exciton, x = e, 
the electron, and x = h, the hole) from the site i to its nearest-neighbor occupied site j is 
defined as 

1% = ^- when Ej-Ei<0 (9) 

^~hop 

and 

1% = 4- expi-fxAEij/kBT) when E j - E { > (10) 

^~hop 

where AEij = Ej — E^. When the site j is not occupied, as it may occur in two- and 
three-dimensional systems, then Pf, = 0. Here, < f x < 1 is an additional parameter 
related to the band structure of nano-crystals, and t£ is a characteristic time for tunneling 
of x— carriers between nearby (n.n.) Si nano-crystals, which is temperature independent. In 
Eqs. H and [lOj, no dependence on the intersite distance, r^, between n.n. Si nano-crystals is 
explicitly shown. In a more refined version of the present model, one could consider in Eqs. 
|9] and [H] a dependence on of the form 

-j- = v x exp ( ) (11) 

T hop 
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where u x is an attempting frequency and 7^ is related to the energy barrier which separates 
the QDs at site i and j. The quantities 7^ and may then be considered as random 
variables. For simplicity, we use r£, p = r^ op here, independently of the type of carrier. Thus, 
the dependence of hopping rates on x is only contained in the factor f x . 

In our calculations, we have arbitrarily chosen f e = 1/3, fh = 2/3 and f ex = 0.4 in order 
to take into account the distribution of the energy-gap discontinuity, E — Eeuik, among 
the valence and conduction bands. Extensive Monte-Carlo simulations have been done 
by using the same set of parameters but for different f ex . The effect of variing f ex is to 
change the barrier for hopping, f ex AEij, which is mostly relevant at low temperatures. 
By increasing f ex , a decrease in the luminescence decay rate occurs, but essentially no 
changes in the luminescence decay lineshapes are observed, i. e. in Eq. p] r increases and j3 
remains constant. The mean distances over which the excitons move before recombination 
are reduced as the effective barrier, f ex AEij, is increased. However, for f ex in the interval 
[0.2 — 0.6] the numerical results are qualitatively not affected and quantitatively they change 
by less than 10 %. 

Before going on with the discussion of the results, we briefly summarize the free param- 
eters entering the model: 

(1) E cut , absolute energy difference cut-off of nearby occupied sites in the lattice. 

(2) Tg^, prefactor for the singlet radiative recombination time (Eq. ||). 

(3) t$, prefactor for the non-radiative recombination time (Eq. ||). 

(4) Thop, prefactor for the tunneling of carriers between nearest-neighbor occupied sites (Eq. 
|and 0). 

Table I reports the values of all the parameters used in our simulations for ID, 2D and 
3D systems whose results are presented in the following. 
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D. The case of isolated Si nano-crystals 



Let us consider the case in which no hopping between n.n. Si nano-crystals occurs in the 
system. This limiting situation may well describe p-Si samples in which Si nano-crystals are 
efficiently passivated and surrounded by thick surface oxide layers. This case represents an 
ensemble of isolated QDs, each characterized by its own size d and transition energy E. To 
obtain the luminescence spectrum of the whole ensemble as a function of time, for different 
observation energies and temperatures, the rate equation describing the time evolution of 
the density N(t) for a single carrier at a given Si nano-crystal (site) is solved. Such a rate 
equation reads, 

dN(t) N N 



dt Trad Ti 

and the exact solution is 



(12) 



N(t) = N(0) exp(-t/r) with - = — + —. (13) 

Typical behaviors of the different times entering in Eq. [13] as a function of temperature, for 
fixed energy E, are shown in Fig. |2]. 

The luminescence L(t) from a single Si nano-crystal is just L(t) = A^(t)/r ra( j, i.e. 

L{t) = L(0) exp(-f/r) (14) 

yielding simple exponential decay, i.e. (3 = 1 (cfr. Eq. [[]). As pointed out in Ref. |§, 
the luminescence decay time is essentially determined by non-radiative recombinations at 
room temperature. According to Eq. [13] and Fig. [| non-radiative recombinations play a 
main role also at low temperatures when the radiative lifetime, dominated by the triplet 
recombination lifetime, becomes very large. In the intermediate temperature range, when 
the quantum efficiency also increases, the singlet lifetime may dominate the recombination 
lifetime, if it becomes sufficiently small and thermal quenching of the phonon population 
takes place, decreasing the capture cross-section to the non-radiative recombination centers. 
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This is case 3 in Fig. |2|. At sufficiently high temperatures, non-radiative recombinations 
dominate again. 

The total radiated intensity I(E,T) = j °° dt L{t) = N(0) r/r TSL d, is related to the 
corresponding quantum efficiency, e(E,T) defined as 

e(E,T) = I -^ = —^— . (15) 
(0) r nr + r rad 

For illustration, the quantity e(E,T)P(E), representing the quantum efficiency of the en- 
semble of QDs described by Eq. [|for an energy E, is plotted for several temperatures versus 
energy in Fig. [3[ The total quantum efficiency of the ensemble, e(T), is just the sum of Eq. 



15] over the total number of Si nano-crystals in the system, i.e. e(T) = J2e P{E) e (EjT). 
This quantity is shown in Fig. |J At high temperatures, we find an exponential behavior 
e(T) ~ exp(— T/Tq), where T = 59 K (see cases (1) and (2) in Fig. in good agreement 
with recent experimental results (see e.g. Ref. pifl ). The value Tiu) p h/kB = 800 K and 
p = 25 (see Section II. B), were chosen such that e(T) displays a maximum in the range 
100 K< T < 200 K and a value of about 10 % at room temperature, as is experimentally 
observed. |5J 

Finally, the absorption coefficient a(E) of the ensemble can be obtained from the relation 

a(E,T)^ C -^\Mf P (E) } (16) 

where \M\ is the dipole matrix element at the band-edge describing radiative transitions 
in the QD and p(E) is the total density of states. \M\ is related to the oscillator strength 
of the transition and hence to the radiative lifetime, i.e. |M| 2 ~ t^ l (E,T). p8| Using 
the result obtained in Section II. A for p(E), Eq. f|, and the expression Eq. || for r ra d, we 
obtain the absorption coefficient of the ensemble as displayed, for various temperatures, in 
Fig. [5]. Remarkably, this very simple model is able to reproduce the main features of the 
measured absorption coefficient, i.e. a quite sharp edge at about 1.6 eV and an increasing 
exponential shape at high energies with a slope which depends on the temperature. An 
exponential fit at high energies, a(E) ~ exp(E / Eq) , of the curves shown in the figure yields 
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for the energy slopes: E = 0.55 eV for T = 30 K, 0.34 eV for T = 100 K, 0.29 eV for 
T = 200 K and 0.28 eV for T = 300 K. The temperature dependence of a(E) is due to 
the r ra d term in Eq. |16[ The exponential dependence of the absorption coefficient has been 
repetitively observed in p-Si [^,0 an< ^ attributed to an Urbach-like edge as in amorphous 
Silicon. The origin of the Urbach-tail is not precisely known so far but it is believed to 
be due to disorder. In the present model, it results from the overlap of the 0-dimensional 
densities of states associated to the QDs randomly arranged in space. 



E. Rules for the Monte-Carlo simulations 

In the case in which hopping of carriers between nearby Si nano-crystals takes place, the 
rate equation for N(t), Eq. [12|, is modified to 

— it 1 = --r--- L -H Pa N i(t) + E Pa N At) (17) 

'rad ' nr j j 

yielding a system of coupled differential equations for the different occupied sites i. For 
simplicity, we have omitted the index x, denoting the type of carrier, from Eq. [T7|. The 
third term in Eq. |17] represents the outgoing particles from site % (loss term), while the 
last term the incoming particles to site i from the neighbor sites j (feeding term). In one 
dimension, the number of nearest neighbor sites n — 2, and because of structural disorder, 
in two dimensions 1 < n < 4, and in three dimensions 1 < n < 6. The system in Eq. |17] is 
conveniently solved by Monte-Carlo (MC) simulations when the total number of effectively 
coupled occupied sites % becomes large. In particular, certain constraints for the occurrence 
of recombination processes, such as the simultaneous presence of an electron and a hole at 
a given site, can be easily implemented with the MC method. 

Before discussing the MC-rules, we need to determine the unit of time, denoted as To, 
which should be sufficiently small such that faster transition events are well described. Once 
the transition times r have all been determined in the system, one can take To = Tf ast /nt, 
where Tf ast is the smallest transition time in the ensemble and n t > 1. We have used n t = 10 
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in our simulations. Notice that the time t becomes now a discrete variable, i.e. t = nr , 
where n > 1 denotes the nth MC-step. The total elapsed time for each MC-step is just r . 

We can now discuss the MC-rules for the present model. Since we are interested only in 
the case of very low carrier density i.e. different electron-hole pairs in the system do not 
see each other (no interaction effects or Auger recombinations), we study the time evolu- 
tion (trajectory) of a single electron-hole pair (either as an exciton or as two independent 
particles), which is initially located (t = 0) at the center of the lattice. Averages are then 
performed over many trajectories and different realizations of disorder. 

During the nth MC-step, a particle at site i may undergo one of the following four 
different processes: 

(1) decaying radiatively (annihilate) with a probability p Ta d = W r rad? 

(2) decaying non-radiatively (annihilate) with a probability p nr = r /r nr , 

(3) hopping to a n.n. site j with a probability = r Pij, or 

(4) remaining at site % with a probability p sit = 1 — p Ta ,d — TV — J2j Pij- 

In the case in which an electron and a hole are considered, a decay can only take place 
when both carriers are located at the same site. 

According to our definition of r , all the above probabilities are smaller than one as 
required. To decide which event will take place for a given particle at site i, we generate 
a random number r, uniformly distributed in the range < r < 1. Then, we evaluate the 
partial sums over the probabilities, i.e. J2T=i Vk = S m , where the index k = 1 represents a 
radiative recombination event, k = 2 a non-radiative one, etc., and m > 1. The successful 
event is the one for which S m > r for the first time. For example, consider a ID system 
and let p ra d = p nT = 0.2 (equally likely radiative and non-radiative recombinations) and 
Pi(j-i) = Pi(i+i) = 0.1 (i.e. equal hopping rates to two n.n. sites). Then, if r = 0.55, 
the particle will undergo the fourth possibility, i.e. it will hop to the second n.n. site. If, 
however, r > 0.6, the particle will remain sit at its present site. 

A single MC-step is completed when the above discussed procedure has been applied 
to each particle in the system (either the exciton, or the electron and hole), and the time 
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t is increased by To, t = t + To. At the nth MC-step, corresponding to time t = nro, the 
luminescence of the system is simply obtained by counting the total number of radiative 
events during the nth MC-step. Since we are also interested in the energy dependence of the 
luminescence, those radiative events originated at sites characterized by transition energies 
E' close to the observation energy E (i.e. \E' — E\ > 5E) are recorded separately. In our 
simulations we used 8E = 0.01 eV. 

We proceed with the numerical results obtained for one-dimensional lattices. 



III. RESULTS FOR ONE-DIMENSIONAL LATTICES 

At very high porosities, p-Si is characterized by a filamentary structure having a quasi- 
one dimensional character on small length scales. |33j Thus, we start our discussion about 



transport (diffusion) of carriers in p-Si on a simplified one-dimensional model. Preliminary 
results with the set of parameters named case (3) in Figs. |2] and ^ are reported in Ref. ||19| . 



When the carriers are allowed to hop between n.n. sites of the linear lattice, the diffusive 
motion takes place in a highly irregular energy landscape. This feature of the model is 
illustrated in Fig. where the local transition energies E are plotted as a function of 
position on the one-dimensional lattice. In the figure, three cases are considered for selected 
values of the parameter E CVLt (see Section II. A). Although local variations in energy can 
be smoothed out when E cut becomes sufficiently small, on large length scales appreciable 
variations of E will always occur (see Fig. |6|). This has important consequences when 
considering hopping processes, because carriers must eventually overcome an effective large 
barrier to diffuse out of a region of local minimum energy. According to our definition 
of the transition rates in Eq. [10], thermally activated hopping against such large scales 
effective barriers will be strongly hindered at sufficiently low temperatures, even in the case 
of low E cut values. The low energy sites with low p ra( j and high energy barriers AEij act as 
temporary traps. Excitons which relax to these sites are temporary trapped into them and 
have to wait long times before being released. 
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To quantify the effect of disorder on diffusion we calculate the mean square displacement 
of the carrier as a function of time, R 2 = (r 2 (t)) — (f(t)) 2 , and investigate its temperature 
dependence. f(t) is the displacement with respect to the initial position, which is assumed 
to be the origin. To study the effect of the dynamics of carriers on the time dependence of 
L(t), we have solved Eq. [T7] for the values given in the Table I. Since in experiments all sites 
of the system are homogeneously excited at t = 0, including those with energy E close to 
the observation energy -E bs, initially, the electron- hole pair is placed at a site (the center of 
the lattice) characterized by an energy E = -B bs- 

In the case of excitons, the behavior of the luminescence is shown in Fig. |7| for four dif- 
ferent temperatures and observation energy E^ = 1.86 eV. Two remarks are evident: i) the 
decay departs from the simple exponential decay obtained for isolated QDs (corresponding 
to Tho P = oo), and ii) the time scale of the decay is strongly temperature dependent. The 
last item is due to the thermal population of the singlet state (fast radiative recombination 
time) and to the onset of efficient non-radiative decay, especially at the highest temperature. 
The lineshape analysis shown in Fig. evidences that a stretched exponential function is 
required to fit the decay. The single exponential or the double exponential functions are 
able to fit the decay only on a limited time range. 

The stretched exponential decay reflects the fact that a time delay between recombination 
events at a given QD may occur when the exciton can diffuse out of the QD. Such time delay 
can persist on large time scales since the exciton may either diffuse far away from the site 
at which the radiative event is expected, or may become temporarily trapped at a nearby 
site with a lower transition energy. When the observed radiative events are collected from 
sites with similar transition energies E, the random local environment around each of these 
sites leads to a distribution of recombination times, which are in general different than the 
single value t T3 ^(E, T) for that energy E. This dispersion of recombination times causes the 
stretched exponential behavior observed in the model. 

In the case of uncorrelated electron and hole motion, a quite different behavior is observed 
(see Fig. ||]). At high temperatures, and intermediate times, an approximate power-law decay 
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of the form 



L(t) ~ t~ a (18) 

takes place. The range of times over which the power-law decay seems to occur grows as 
the temperature is raised. Here, the uncorrelated motion of the electron and hole leads to 
a dramatic slowing down of the luminescence decay, since the two carriers are no longer 
constrained to be all the time simultaneously at the same site. Similar lineshapes have 



been observed for amorphous hydrogenated Silicon. |34| In this system, the luminescence is 



explained as due to geminate recombination between e-h pair localized in band tail states. 



We have performed extensive MC simulations for different temperatures. In the case 
of excitons, the decay lineshapes are fitted by a stretched exponential function for all the 
temperatures considered. A statistical weigth of the simulation points has been used to 
obtain the best fit function. The results of the fits are summarized in the upper panel of 
Fig. |TU[ Typical errors for r of 10 % and (3 of 0.02 are expected, due to statistical fluctuations 
and goodness of the fits. The lifetime r decreases significantly as the temperature is raised, 
while the dispersion exponent /3 increases and tends to unity at high temperatures. 

We have also calculated the quantum efficiency e(E,T) as a function of E for a fixed 
temperature, displaying a maximum around E = 1.86 eV and found that hopping of carriers 
does not affect sensitively e(E, T). This is clear because e(E, T) is essentially determined by 
the on-site recombination dynamics. Independently of the dispersive motion, all the carriers 
recombine radiatively or non-radiatively after a sufficiently long time. The diffusion only 
influences the recombination dynamics by introducing long lived recombinations. 

In Fig. |11|, we show the values of R 2 (t) for different temperatures. At intermediate and 
high temperatures, an anomalous behavior of the form R 2 (t) ~ t 2 ^ dw with d w > 2, takes 
place. Such anomalies in R 2 (t) are typical of diffusion on fractals, where the presence of 
dangling ends and loops on all length scales in the fractal slows down the diffusion process on 
all time scales. On fractals, however, the diffusion exponent d w is temperature independent, 
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3~6fl while here d w turns out to depend strongly on temperature (see Fig. [T2| ). The anomalies 
in R 2 are due to the random distribution of hopping rates, which is known to yield values of 



d w > 2 and dependent on temperature (see e.g. Ref. |37j|). At low temperatures, however, 
the time dependence of R 2 (t) departs from a power law. 

At sufficiently low temperatures, indeed, the present model is expected to display trans- 
port behavior similar to diffusion in the presence of random fields (Sinai model), see Refs. 
|38| and |39| . To show this, we have drawn in Fig. [13] a small section of a typical energy 



landscape versus position. The arrows represent the local fields felt by the carrier which 
is proportional to the energy difference AEij between the two n.n. sites % and j. Due to 
the assumed random distribution of nano-crystal sizes in the system, the local fields are 
randomly oriented, similarly as in the Sinai model. Diffusion of carriers in the system is 
biased by these random fields and becomes ultra-anomalously slow. Of course, in the present 
model the energy barriers (and also the potential valleys) can not grow indefinitely as in the 
Sinai case. However, when the temperature is sufficiently low, diffusion is strongly hindered 
and the carriers can only explore small length scales since the effect of the random fields 
becomes dominant. In these circumstances, a logarithmic time dependence is expected 

R 2 (t) ~ (logt) 4 (19) 



at long times. Results obtained for the case of excitons are plotted in Fig. [TJ. These 
suggest that for large times and low temperatures, our model displays logarithmic time 
dependencies of the form predicted by Eq. pl| To our knowledge, this is the first time that 
such logarithmic time dependencies of R 2 (t) are predicted for a model aimed to describe a 
disordered system such as porous Silicon. Actually, the value of the exponent describing the 
logarithmic time dependence turns out to depend slightly on temperature, and the value 
of four is observed only for T = 30 K. At lower T, R tends to extrapolate to a constant 
value because non-radiative annihilation events become dominant, and the carriers can only 
explore a finite spatial extent asymptotically. A similar trend is also observed for electrons 
and holes in the case of uncorrelated electron-hole pairs. 



IV. TWO- AND THREE-DIMENSIONAL MODELS AND RESULTS 



The actual geometrical structure of porous Silicon is not known accurately. For many 
purposes, however, one may hope to capture the essential features of such complex structures 



by modeling them with simple, yet non trivial, percolation-like clusters. [36] The clusters 
are generated by using a modified version of the well-known growth algorithm employed for 
percolation clusters (see e.g. Ref. which is adapted here to our present purposes. 

For simplicity we consider a square lattice in two dimensions, and a simple cubic one in 
three dimensions. The linear size of the lattice is denoted by t. The growing process starts at 
the seed, which is located at the center of the otherwise empty lattice. The growth proceeds 
according to the following rules. A nearest neighbor site j of the seed can be occupied 
with probability p n , with the index n indicating the number of occupied sites nearest to the 
growing site j. Initially, n — 1. If the site j is not occupied, i.e. if it does not become part 
of the cluster, it is blocked and cannot be occupied later. The process continues from the 
last occupied sites and now values n > 1 may occur. 

For standard percolation clusters, one takes p n = p independently of n. When p < p CT a, 
p CTit = 0.593 (two dimensions), and p CT - lt = 0.312 (three dimensions), only finite clusters 
can grow, while infinite clusters develop when p > p crit . When p = p crit , large percolation 
clusters may grow which are fractal (with fractal dimension df = 1.896 in two dimensions, 
and df = 2.5 in three dimensions) on large length scales. Finite clusters and the infinite 
cluster above p cr \ t are also fractal (with the same fractal dimension df) for length scales 
smaller than the correlation length. The latter diverges at p = p C rit- 

In our model, we take p\ > P2 and p n = when n > 2. In the following, we consider 
the cases, p\ = 0.65 and P2 = 0.15 in two dimensions, and p\ = 0.45 and P2 = 0.17 in three 
dimensions. 

We notice that when p\ — > 1, and P2 — > 0, the clusters tend to grow linearly, while for 
relatively large P2, compact clusters can be grown. By varying both pi and P2, a variety of 
structures can be obtained which are suitable for our present purposes. A measure of the 
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porosity of the cluster is simply given by the ratio between empty sites and occupied sites. 
Typical examples are shown in Fig. [15] for three dimensional clusters, and in Ref. ^0] for 
two-dimensional clusters. The clusters are uniform on large length scales, but still display 



spatial fluctuations on small length scales as the infinite cluster above p C nt- |20] Thus, the 
intrinsic fractal character of percolation is common to the clusters shown in Fig. O too. It 
is possible to construct highly constrained and filamentary clusters (like that in Fig. |l5|(a)) 
or more compact clusters (like that in Fig. |l5"|(b)). 

We have studied the dynamical behavior of carriers on these clusters, following the 
same MC-rules described above. The question now is how the extra degrees of freedom in 
space, resulting from the two- or three-dimensional topology of the clusters, modify the time 
decay of L(t). Preliminary results for two-dimensional systems have been reported in Ref. 
|40| . Results for three-dimensional systems are reported in Fig. [16] for the luminescence 
decay. pOJ For excitons and uncorrelated electron-hole pairs, the lineshape is a stretched 
exponential and a power law, respectively. Hence the decay lineshape is not modified by the 
increased dimensionality. Results of stretched exponential fits to the luminescence decay in 
two or three dimensions, for the exciton model, are reported in Fig. [10 as a function of 
temperature. 

The three-dimensional topology has important consequences on the decay of the lumi- 
nescence. The t and (3 values are lower, indicating that the role of hopping processes in 
three dimensions is more important than in lower dimensions. While the one-dimensional 
model (and also the two-dimensional one) predicts values of (5 close to unity already at room 
temperature, in three dimensions, in contrast, the theoretical values are consistent with the 
experimental results for (3, which typically saturate at values (3 = 0.7 — 0.8. |18],|19j Thus, 
values of /3 < 1 at high temperatures, and intermediate porosities, can be explained by our 
model as a result of the interplay between a complex conducting matrix (representing the 
topology of p-Si), and an additional local disorder due to the distribution of nano-crystal 
sizes in the system. 

A comparison of the temperature dependence of r and j3 for two different sets of pa- 
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rameters (but for the same growth parameters pi, P2 mentioned above) is reported in Fig. 



IT , The effect of variing r nr is evident in (3: the larger the r nr the lower is (3. The role of 
temporary traps, which tend to reduce (3 and are effective at low temperatures, is enhanced 
when r nr is increased due to a relative increase of the role of exciton diffusion (~ r nr /rho P ). 
The dynamics of carriers is influenced by the parameters p\ and P2 used to construct the 



clusters. The more the cluster is filamentary (see e.g. case (a) in Fig. |15|), the more the 
excitons are constrained to move along one direction, and the larger are r and (3. Then, the 
dynamical behavior of carriers should have essentially a one- dimensional character, at least 
for sufficiently low temperatures where transport along the structure is strongly hindered. 
Thus, we expect to observe, in this regime, the same logarithmic time dependence of R 2 {t) 
as in the one-dimensional model, Eq. [19[ at short times. 

For arbitrary porosities and not too low temperatures, we find anomalous behavior of 
R 2 (t) ~ t 2 l dw , at intermediate times (Fig. pl|). Here, d w > 2 and increases as the tempera- 
ture decreases, more than in the one-dimensional case (Fig. |I2"D . Eventually, a logarithmic 
time dependence of R 2 (t) is observed for sufficiently low T. In three dimensions, we have 
found clear evidence for the behavior 

R 2 ~(\ogt) 2 (20) 

at long times (Fig. [19]), in contrast to the fourth-power result valid for the Sinai model. This 
new behavior represented by Eq. [20 has been predicted for the case of biased diffusion on 



random fractals in two dimensions. [39| As we can see, a similar behavior can be expected 



also in three dimensions. This feature remains to be studied theoretically. 



V. DISCUSSION 

According to the present model, the photoluminescence decay, L(t), is determined, in 
the case of excitons, by the recombination and the escape rates from the isolated QD. In the 
case of uncorrelated electron-hole pairs, L(t) is determined by the diffusion of the electron 
and the hole, because both need to be on the same QD to recombine radiatively. 
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In the case of excitons, the disordered environments of the QDs, from which the lumi- 
nescence is observed, cause a distribution of waiting times for hopping and/or release times 
for activated emission from the temporary traps, and the luminescence displays a stretched 
exponential decay. In the case of uncorrelated electron-hole pairs, the main role is played 
by the encounter probability of the electron and the hole to be on the same site. This is 
determined by the diffusion of the slower moving particle, i.e. the hole, and the time de- 
pendence of the recombination rate becomes a power-law, which is then reflected in the 
luminescence decay. 

For all spatial dimensions D, the increase in (3 is almost linear for 10 <T<60 K, has 
a plateau for 60 <T < 200 K and then increases rapidly to 1 at higher temperatures. 
This behavior corresponds to the following three typical regimes: a first one in which the 
temporary traps play a role (low temperatures), a second one in which hopping dominates 
and diffusion is restricted only by the p-Si network geometry (intermediate temperatures), 
and a third one, in which the rapid on-site non-radiative recombinations are most effective 
and dominant (high temperatures). 

Whereas there is general agreement among different experimentalists that r strongly de- 
pends on temperature, |7|,|,[!1| the situation is less clear regarding the exponent (3. While 



some authors report values of (3 = 0.8 — 0.9, independent of the temperature |15|,[16l , others 



found temperature dependent values. |17Hl9|l In this case, the (3 values follow similar tem- 
perature and energy dependences in the range 0.4-0.8. By varying the parameters in our 
model, we are able to describe both situations. Temperature dependent (3 values are shown 
in Fig. |17[ Temperature independent values can be obtained when Th op < T ni andr S i ng , since 
in this case the activation term in Eq. [10], which is mainly responsable for the temperature 
dependence of (3, does not influence the diffusion of carriers, being only restricted by the 
geometrical constraints of the p-Si network. 

By going from an interconnected array of QDs to an ensemble of isolated QDs our model 
predicts that the lifetime r and the dispersion exponent (3 increase, with (3 taking values near 
1. This has been indeed observed in several experiments. For example, the data presented in 
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Ref. [p] , and reanalised using stretched exponential functions, show that r and (3 increase 



after dry oxidation of p-Si. It is well know that the interconnected array of QDs typical of 
p-Si is transformed after dry oxidation into a dispersion of nano-crystals immersed in a SiO x 
matrix; increasing further the oxidation a porous glass is formed. U 

Our simulations show that hopping of excitons is responsable for the stretched exponen- 
tial decay observed in p-Si. This rules out other explanations which are based on analogies 
with hydrogenated amorphous Si (a-Si:H), see e.g. Ref. [^2| . In fact we demonstrated that 
PL decay lineshapes, strongly resembling those measured in a-Si:H, are found for uncorre- 
lated e-h pairs. This behavior can possibly explain the lack of room temperature lumines- 
cence in a-Si:H. In fact, it is possible that the fast diffusing species in the pair is rapidly 
trapped into non-radiative recombination centers when the thermal energy is sufficiently 
large to promote the particle to extended states. J43| 



Annealing treatments induce a quenching of the room temperature luminescence with a 
decrease in the values of r and (3. |T3| Within our model, this is simulated by diminishing 
t$ and Thop yielding lower values of r, (3 and quantum efficiency (see Table II). The separate 
effects of variing t$ and r hop are the following: 

(a) Reducing the non-radiative lifetime, yields an increase of j3 and a decrease of r. The rea- 
son is that the re-population of the target QDs, those from which the luminescence decay is 
recorded, is reduced by the competitive non-radiative recombination channels. Less excitons 
reach the target QD at long times because a large part recombine non-radiatively as they 
diffuse through the system. The target QD effectively behaves as being more isolated (in- 
crease of (3). Then, r is reduced due to the competing on-site non-radiative recombinations, 
and consequently the quantum efficiency is also reduced. 

(b) Reducing T\ mp yields a reduction of r and (3. In fact, for higher hopping probabilities, 
the target QDs are emptied faster due to the competing loss of excitons caused by hopping 
(reducing r), while at long times more excitons reach the target QD giving rise to long- 
lived recombinations (reducing (3). In this case, no variations in the quantum efficiency are 
expected. 
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Hence to obtain the trends measured during the annealing experiments both the lifetimes, 
t$ and r hop , should be reduced. 

Let us compare our simulations with other models. The non-exponential decay of L(t) 
has been explained as due to a distribution of r ra d values arising from a shape distribution 
of quantum dots with the same emission energy. ]7J Within this assumption, however, it 
is difficult to explain the temperature dependence of (3 already discussed. It has also been 
proposed that the non-exponential behavior of L(t) is a consequence of a distribution of non- 
radiative decay rates. |44"l , f45f This fact is not considered in our model, and a refined version 
of it should certainly contain such a dependence. However, we want to emphasize here that 
it is the dispersive motion of excitons which is mostly responsible for the non-exponential 
decay of L(t). The radiative and non-radiative lifetime distributions could indeed concur 
to produce a similar effect on L(t), but in isolated and well passivated quantum dots it is 
reasonable to expect essentially a simple exponential decay of L(t). Such a behavior has 
been indeed observed for some p-Si samples. |56| 

The strong temperature dependence of the effective anomalous diffusion exponent d w , 
and the logarithmic time dependencies of R 2 {t) at lower temperatures, have important conse- 
quences on the ac-conductivity, which is expected to display a strong frequency dependence 
at low temperatures. This observation is supported by recent experimental results on the 



ac-conductivity of p-Si. |47| In contrast, at higher temperatures, diffusion becomes normal, 
R 2 (t) ~ t, i.e. d w — > 2, and the dispersion frequency range shrinks considerably. It should 
be emphasized that a model for the structure of p-Si has been suggested recently, [[47 



in 



which the p-Si network is assumed to be the infinite percolation above criticality. From 
our present results, however, one can see that such a model is far too simple to describe 
p-Si samples of different porosities, and a more elaborated model is required. In addition, 



theoretical results for the ac-conductivity discussed so far, |TJ] have been obtained using 
mean-field approximations, and the new intrinsic behavior suggested here for diffusion of 
carriers in p-Si at low temperatures can not be obtained. 
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VI. CONCLUSION 



The theoretical results obtained in this paper have for the moment a semi-quantitative 
character. A quantitative comparison between the present results and the available ex- 
perimental data can be made possible when the parameters entering the model may be 
estimated independently and more accurately Also, the precise relations between the nano- 
crystal structure and dynamical properties are so far unkown. We have tried to fill this 
lackness by making ad hoc, yet standard assumptions which seem to describe the physics 
of p-Si rather well. In addition, the actual values of the parameters are so strongly sam- 
ple dependent and influenced by various treatments (ageing, oxidation, storage, excitation 
conditions, etc.) that quantitative predictions valid in general are not possible. However, 
several experimentally measured trends in the photoluminescence lifetimes, dispersion ex- 
ponents and quantum efficiencies of the luminescence are correctly explained in the frame 
of the present calculations by choosing the free parameters of the model appropriately. This 
gives a strong support for the validity of the present model for describing the recombination 
dynamics in p-Si. 

Finally, by assuming different forms for the various quantities reported in Section II, other 
systems can be described as well by the present model, as e.g. nanometer sized Si crystallites, 
P8fl CdSSe quantum dots, jKJ etc. All these systems show a stretched exponential decay of 
the luminescence. 
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TABLES 

TABLE I. Values used in the Monte-Carlo calculations. ft^/D) j g ^ e \[ neav s [ ze Q f the 
D-dimensional lattice employed, t max the maximum number of MC step for a single run, n mns 
the number of runs (i. e. the number of different particles used) for a single cluster realization, 
n CO nf the number of different cluster realizations over which the reported results are averaged. The 
meaning of the other parameters is given in the text. All the results presented in this paper refer 
to this choice of parameters unless otherwise stated. 
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TABLE II. 300 K simulation results for a 3D clusters of 65% porosity as a function of the time 
parameters: Tn? and Th op . The other parameters are given in Table I. 
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FIGURES 

FIG. 1. Distribution of optical-band gaps for an ensemble of quantum dots (dotted line) (Eq. 
2), and corresponding density of states (continuous line) (Eq. 4) versus energy. 

FIG. 2. Temperature dependence of the transition times (Eq. 13) for an energy of 1.86 eV: r ra d 
(short-dashed lines), r nr (long-dashed lines) and r (continuous line). Here we have used the same set 
of parameters as given in Table [VT| except for: case (1), thick lines: = 150 fis, ri? = 900 /js; 
case (2), thin lines: r s ^ g = 100 [is, rffl = 400 [is; and case (3), thin lines: r s ^ g = 200 /is, 
T £ ] = 7000 fjs. 

FIG. 3. Quantum efficiency of an ensemble of isolated quantum dots: e(E, T)P(E) versus E for 
various temperatures reported on the curves in K. The temperature dependence of the Si band-gap 
has been neglected. 

FIG. 4. Total quantum efficiency of the ensemble of quantum dots, e(T), versus temperature. 



Here we have used the same set of parameters as given in Table VI except for: case (1), thick 



line: r s ^* g = 150 /is, t$ = 900 [is; case (2), thin line: r s ^ g = 100 [is, Tn? = 400 [is; and case (3), 
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thin dashed line: r ^„ = 200 /is, t$ = 7000 /is. The thick dashed lines represent the exponential 
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behavior e(T) ~ exp(— T/Tq), with To = 59 K for case (1) and case (2) and To = 75 K for case (3). 

FIG. 5. Absorption coefficient of the ensemble of isolated quantum dots, a(E, T), versus energy 
for the indicated temperatures, T. 

FIG. 6. Snapshot of the energy landscape for diffusion in a one-dimensional lattice. The tran- 
sition energy E is plotted versus position, for: E cut = oo (uncorrelated values, top panel), 0.04 eV 
(middle panel), and 0.004 eV (bottom panel). The mean energy value over the 29791 sites is 1.86 
eV for the three cases. 

FIG. 7. Monte-Carlo simulation results for a one-dimensional system of excitons. Luminescence 
versus time for four different temperatures: T = 20 K (squares), 100 K (triangles), 300 K (discs) 
and 500 K (empty circles). 
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FIG. 8. Fitting of the luminescence decay (circles) obtained from a Monte-Carlo simulation for 
a one-dimensional system of excitons at T = 100 K. The following fitting functions and best fit 
parameters have been used: Stretched exponential function with r = 345 /jls and (3 = 0.815 (full 
line); single exponential decay with r = 555 fis (dotted line); double exponential with n = 80 ^s 
and T2 = 503 /xs (dashed-dotted line). The top panel reports the relative errors, defined as the ratio 
of the difference of the simulation data (Y) minus the computed value (Yfit) times the simulation 
data, i. e. [(Y — Yfit) /Y], for the three fitting functions. 

FIG. 9. Monte-Carlo simulation results for a one-dimensional system of uncorrelated electron 
and hole pairs. Luminescence versus time for four different temperatures: T = 30 K (dots), 50 K 
(circles), 100 K (triangles) and 300 K (squares). The lines are power-law fits to the decay with the 
following exponents: a = 0.86 (T = 30 K), 0.79 (T = 50 K), 0.73 (T = 100 K) and 0.91 (T = 300 
K). The curves are vertically shifted for clarity. 

FIG. 10. Temperature dependence of the luminescence recombination time r (empty discs) and 
dispersion exponent (3 (filled discs) obtained by a stretched exponential fit to the time decay of 
the luminescence, for an exciton population in one-dimensional (upper panel), two-dimensional 
(middle panel) and three dimensional (lower panel) clusters. 

FIG. 11. Mean square displacement R 2 (t) of excitons in one dimension for three different 
temperatures: T = 20 K (circles), 100 K (triangles) and 500 K (squares). Note the departure from 
a power-law behavior at low temperatures. The curves are vertically shifted for clarity. 

FIG. 12. Anomalous diffusion exponent d w obtained by a power-law fit to the mean 
square displacement R 2 (t) of excitons in one-dimensional (circles), two-dimensional (discs) and 
three-dimensional (triangles) clusters. 
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FIG. 13. Energy landscape (transition energies E vs position) in one dimension. The arrows are 
drawn in direct correspondence with the local energy difference between nearby sites. The resulting 
one-dimensional chain (lower panel) is equivalent to the Sinai model, in which uncorrelated random 
fields are present at each site of the lattice. 

FIG. 14. Mean square displacement R 2 (t) of excitons in one dimension plotted versus (logt) 4 
for the temperatures: T = 20 K, 30 K and 50 K. The straight line through the T = 30 K results 
represents the asymptotic (logi) 4 behavior. The line through the T = 50 K results shows the 
power law t 2 ^ dw behavior, with d w = 2.5. 

FIG. 15. Percolation-like clusters grown on a simple cubic lattice representing porous Silicon 
of the same 97.4 % porosity but different growth probabilities: (a) p\ = 0.9, p2 = 0.05035, (b) 
px = 0.3235, p 2 = 0.4. 

FIG. 16. Luminescence decay of a three-dimensional system of excitons (upper panel) and of 
uncorrelated electron-hole pairs (lower panel) for three different temperatures: T = 300 K (circles), 
200 K (discs) and 100 K (triangles). The lines are fits to the luminescence decays. Stretched 
exponentials have been used for the exciton case with the following parameters: r = 175 /xs and 
(3 = 0.58 (T = 100 K), r = 135 lis and (3 = 0.62 (T = 200 K), and r = 66 lis and (3 = 0.73 
(T = 300 K). Power laws, t~ a , have been used for the uncorrelated electron-hole pairs with the 
following parameters: a = 0.96 (T = 100 K), a = 1.10 (T = 200 K) and a = 1.28 (T = 300 K). In 
this case, the curves are vertically shifted for clarity. 

FIG. 17. Temperature dependence of the luminescence recombination times r (empty symbols) 
and dispersion exponents (3 (filled symbols) for excitons in three dimensional clusters and for two 
different sets of parameters. Circles correspond to case (1) of Fig. || and to the standard set of 
parameters reported in the Table VI. Triangles correspond to case (2) of Fig. |2| and hopping time 



Thop = 100 ijls. The lines are the lifetimes calculated using Eq. |l3| in the case of isolated quantum 
dots, i. e. without hopping, for case (1) (full line) and case (2) (dotted line). 
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FIG. 18. Mean square displacement R 2 (t) of excitons in three dimensional clusters for four 
different temperatures: T = 20 K (squares), T = 30 K (triangles), T = 50 K (circles) and T = 100 
K (discs). The curves are vertically shifted for clarity. The lines are power law fits. Note, however, 
the bending of the data at low temperatures for short times. 

FIG. 19. Mean square displacement R 2 (t) of excitons in three dimensional clusters for the 
temperature T = 30 K, plotted as a function of (log£) 2 . Note the linear increase at long times. 
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